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Abstract. Using DIRBE and FIRAS maps at high lati- 
tude > 10°) we derive the spatial distribution of the 
dust temperature associated with the diffuse cirrus and 
the dense molecular clouds. For a v 2 emissivity law, we 
find that the equilibrium dust temperature of the cirrus 
is about 17.5 K with only small variations over the high 
latitude sky. Comparison of the far Infrared DIRBE maps 
shows the presence of a colder emission component with a 
temperature around 15 K, assuming a v 2 emissivity law. 
The lowest values of the temperature found in the cold 
regions (~ 13 K) are compatible with the results recently 
obtained for dense cores in star forming regions by the 
balloon-borne experiment SPM-PRONAOS (Ristorcelli et 
al., 1996, 1998, Serra et al., 1997). This cold component is 
in particular present in the direction of known molecular 
complexes with low star forming activity such as Taurus. 
The association between the cold component and molec- 
ular clouds is further demonstrated by the fact that all 
sky pixels with significant cold emission have an excess IR 
emission with respect to the high latitude IR/HI correla- 
tion. We have deduced a threshold value of the column 
density, N^/=2.5 10 20 H cm~ 2 , below which cold dust is 
not detected within the FIRAS beam of ~ 7°. We have 
re-examined the problem of the existence of a very cold 
dust component (T ~ 7 K) by combining DIRBE maps 
of the cold emission with FIRAS spectra, corrected for 
the isotropic component found in Puget et al. (1996). The 
warm and cold component deduced from the analysis of 
DIRBE maps account for the Galactic FIRAS spectra with 
no need for a very cold component (T ~ 7 K). 

Key words: ISM: clouds, ISM: dust, extinction, Infrared: 
ISM: continuum, ISM: structure, ISM: general, Radio- 
continuum: ISM 
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1. Introduction 

Much has been learned about the emission of the inter- 
stellar dust from IRAS observations. To account for the 
Galactic energy emitted in the near to far infrared, it is 
necessary to have a broad dust size distribution from large 
grains down to large molecules. For example, Desert et al. 

(1990) (see also Draine & Anderson 1985, Puget et al. 
1985, Weiland et al. 1986, Siebenmorgen & Kriigel 1992 
and Dwek et al. 1997) have proposed a consistent interpre- 
tation of both the infrared emission in diffuse HI clouds 
and the interstellar extinction curve using a model with 
three components: PAHs (Policyclic Aromatic Hydrocar- 
bons), very small grains (VSGs) and large grains. PAHs 
and VSGs are small enough (a < 10 nm) to experience 
significant temperature fluctuations after photon absorp- 
tion. They emit over a wide range of temperatures and 
dominate the emission for A < 60 /im. The large grain 
component is the more traditional dust component his- 
torically inferred from optical studies. These grains are in 
equilibrium with the incident radiation field with a tem- 
perature of about 17 K in the diffuse atomic medium. The 
large spatial variations of the infrared spectrum over the 
wavelenghts range 12 — 100 /im have been interpreted as 
changes in the abundance of small grains (Boulanger et 
al. 1990, Bernard et al. 1993). 

With the DIRBE (Diffuse InfraRed Background Ex- 
periment) and FIRAS (Far InfraRed Absolute Spec- 
trophotometer) instruments on board the COsmic Back- 
ground Explorer (COBE) satellite, we have a measure of 
the whole emission spectrum of the interstellar dust from 
the near infrared to millimeter wavelengths. The mean 
Far Infrared (FIR) and submillimeter (submm) spectrum 
of the whole Galaxy was first derived by Wright et al. 

(1991) using the FIRAS data. They showed that this spec- 
trum can be fitted by a single temperature component 
with an emissivity index (a) equal to 1.65 (T=23.3 K) 
or by a two temperature model assuming a=2 (Ti=20.4 
K and T2=4.77 K), with a statistically better agreement 
for the two component model. More recently, Reach et al. 
(1995) have followed this analysis by performing a spec- 
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tral decomposition of the FIRAS data for 146 bins over 
the whole sky. They have identified a very cold Galactic 
component (4-7 K) correlated with the warm component 
(16-21 K). The correlation strongly suggests a Galactic 
origin for the very cold component. The presence of such 
a component questioned our understanding of the physics 
of the FIR Galactic emitters. Only clouds with a central 
extinction A v higher than a few tens could sufficiently 
attenuate the heating rate for classical interstellar grains 
(graphite or silicate). Such high extinction corresponds to 
very dense molecular clouds which are not ubiquitous in 
the Galaxy. The presence of a very cold emission compo- 
nent outside high extinction regions would thus demon- 
strate the existence of a yet unknown dust component 
(e.g. needle-like or fractal grains, or very large particles) 
and/or a departure from a v 2 power-law emissivity due 
for example to a sub-mm spectral feature. However, the 
Reach et al. (1995) analysis is based on total power FI- 
RAS spectra which could include a non-Galactic emission 
component. Using HI data in addition to FIRAS spectra, 
Pugct et al. (1996) have found, in the residual emission af- 
ter the removal of the HI correlated emission, an isotropic 
component which could be the Cosmic Far Infrared Back- 
gRound (CFIBR) due to distant galaxies. At \b\ > 40°, 
at least one third of the emission at 500 /im comes from 
the isotropic component which should strongly contami- 
nate the analysis of faint regions of Reach et al. (1995). 
However, this isotropic component cannot explain the ob- 
served correlation between the warm and very cold optical 
depths. 

Boulanger et al. (1996) and Dwek et al. (1997) have 
derived the emission spectrum of dust using the spatial 
correlation of the FIRAS data with Galactic templates. 
This spectrum is insensitive to any isotropic terms (as the 
CFIBR). The spectrum of Dwek et al. was obtained by 
deriving the slope of the /(A) vs 7(100 /im) correlation di- 
agrams at high latitude. Boulanger et al. have used the 
correlation between the FIR and HI 21cm emission to es- 
timate the spectrum of dust associated with HI gas. Both 
spectra are very well-fitted by a single modified Planck 
curve with a v 2 emissivity law and a temperature of 17.5 
K. No significant very cold component is detected in the 
diffuse medium at high latitude. 

In this paper, we present a statistical analysis of the 
dust temperature in the nearby interstellar medium seen 
at |b| > 10°, using DIRBE and FIRAS data. We show 
that the dust in molecular clouds is colder than in dif- 
fuse atomic gas and argue that the detection of very cold 
dust is due to the combined emission of the isotropic back- 
ground and cold molecular clouds. In Sect. 2, we present 
the data we have used. The map of the large grain tem- 
perature is derived and analysed in Sect. 3. The DIRBE 
data at A > 60/im allow us to build the maps of the ex- 
cess of the FIR/60 /im correlation (Sect. 4a). These excess 
maps trace the cold molecular clouds (Sect. 4b-c) and thus 
correspond to the cold component of the dust emission. 



The cold maps, combined with FIRAS spectra, point out 
the significant difference in the dust temperatures in the 
atomic and molecular parts of the ISM in the whole sky 
(Sect. 5.1). The FIRAS residual emission is discussed in 
Sect. 5.2. A general discussion is presented in Sect. 6. 

2. Data presentation and preparation 

DIRBE is a photometer with ten bands covering the range 
from 1.25 to 240 /im with 40 arcmin resolution (Silverberg 
et al. 1993). We restrict our analysis to A > 60 /im because 
our study is focused on the large grain emission of the 
ISM. We use annual averaged maps because they have 
a higher signal to noise ratio than maps interpolated at 
the solar elongation of 90° (see the DIRBE explanatory 
supplement). 

The FIRAS instrument is a polarising Michclson inter- 
ferometer with 7° resolution and two separate bands which 
have a fixed spectral resolution of 0.57 cm -1 (Fixsen et al. 
1994). The low frequency band (2.2 to 20 cm -1 ) was de- 
signed to study the CMB and the high frequency band (20 
to 96 cm -1 ) to measure the dust emission spectrum in the 
Galaxy. We use the so-called LLSS (Left Low Short Slow) 
and RHSS (Right High Short Slow) data which cover the 
low and high frequency bands respectively (see the FIRAS 
explanatory supplement). 

Since in our study we combine the FIRAS and DIRBE 
data at 140 and 240 /im, we have first to check their con- 
sistency. For that, we convolve the FIRAS spectra inside 
the two DIRBE bands at 140 and 240 /im and the DIRBE 
data (at 140 and 240 /im) with the FIRAS Point Spread 
Function (PSF). The PSF is not precisely known for all 
wavelengths, so we use the approximation suggested by 
Mather (private communication) of a 7° diameter circle 
convolved with a line of 3° length perpendicular to the 
ecliptic plane (Mather et al. 1986). The pixel to pixel 
correlation shows that the difference of responsivity be- 
tween the two instruments is ~0.1% at 140 /im and ~1% 
at 240 /im. 

Before studying the Galactic emission, we have subtracted 
the CMB and its dipole emission from the FIRAS data 
using the parameters given by Mather et al. (1994) and 
Fixsen et al. (1994). To remove the interplanetary dust 
emission, we consider the 25 /im map as a spatial tem- 
plate for the interplanetary dust emission. For the 60 and 
100 /im maps, we use the zodiacal emission ratios given by 
Boulanger et al. (1996). At longer wavelength, we follow 
Reach et al. (1995) by taking a zodiacal spectrum I v oc v 3 . 
To analyse the Galactic emission, it is necessary to choose 
a reference zero emission level because the sky brightness 
may contain an isotropic component (Puget et al. 1996). 
Following Puget et al., we choose as a reference the ex- 
trapolated value for zero HI column density of the FIR- 
HI correlation. We found an extrapolated emission of 0.65 
MJy/sr and 0.56 MJy/sr for the DIRBE 240 and 140 /im 
data (convolved with the FIRAS PSF). These values are 
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lower than the minimum emission measured in the maps 
(equal to 1.70 ± 0.07 and 1.67 ± 0.1 MJy/sr at 240 and 
140 /im respectively), which shows that all DIRBE pixels 
(at 7° resolution) contain a significant Galactic contribu- 
tion. For the analysis of the large grain temperatures, in 
Sect. 3, we subtract from the DIRBE data at 240 and 
140 ^im 0.65 and 0.56 MJy/sr respectively. From the FI- 
RAS data, we subtract the zero HI column density spec- 
trum of Puget et al. (1996) smoothed to the resolution 
of 5.7 cm -1 . It is important to check that the reference 
flux is consistent between DIRBE and FIRAS. For this, 
we compare the DIRBE and FIRAS data at high Galactic 
latitude at 240 /^m and find an offset I(DIRBE)-I(FIRAS) 
equal to -0.15 MJy/sr. The DIRBE reference zero emis- 
sion level at 240 /im corrected for this difference is in per- 
fect agreement with the residual spectrum of Puget et al. 
(1996). The brightness values subtracted from the DIRBE 
data should not be used as an estimate of the CFIBR since 
it may contain some Galactic emission that is associated 
with the diffuse ionised gas. These brightness values are 
lower than the upper limits on the CFIBR (1.76 and 2.66 
MJy/sr at 240 and 140 /im respectively) derived by Hauser 
(1995). 

3. Large grain temperature maps 

In this section, to study the spatial variations of the tem- 
perature, we assume that each individual spectrum can be 
represented as a single modified Planck curve with a v 2 
emissivity law. Obviously the hypothesis of a single tem- 
perature is not valid in molecular regions and near the 
Galactic plane where several clouds with different tem- 
peratures overlap. The temperatures of large grains are 
computed from the ratio of DIRBE emissions at 240 and 
140 /im (taking into account the bandpass) at the FIRAS 
resolution. The very cold component (T ~ 5 K) if it ex- 
ists does not significantly contribute to these bands. We 
have not used the DIRBE data at shorter wavelengths (60 
and 100 /im) since (1) small grains contribute significantly 
to the emission (Desert et al. 1990, Sodroski et al. 1994, 
Laureijs et al., 1996) and (2) the reference zero emission 
level is difficult to evaluate because of the residual zodia- 
cal emission. 

Fig. 1 shows the temperature map. Uncertainties are 
about 5 mK for \b\ < 10°, 50 mK for 10° < \b\ < 30° 
and 300 mK at higher Galactic latitudes. These errors 
are estimated using uncertainties in the DIRBE 7° flux 
obtained by convolving the rms 40' DIRBE maps with 
the FIRAS PSF. 

Most of the variations visible on Fig. 1 are real. Star 
forming molecular regions (e.g. Orion and Ophiucus) have 
a higher temperature due to dust heated by embedded 
young stars, while the temperature of clouds forming low 
mass stars is lower (Taurus, part of Cygnus, Chamaeleon, 
Polaris...). Cold spots in the south and north polar caps 



are not significant since they correspond to regions with 
very low emission in the 140 and 240 /im maps. Cuts along 
the Galactic plane confirm the decrease of the temperature 
with increasing galactocentric distance already studied by 
Sodroski et al. (1994). 

The scatter diagram presented in Fig. 2 does not show 
any systematic variation of the temperatures with the 
Galactic latitude, which suggests a relatively constant In- 
terstellar Radiation Field (ISRF) on large angular scales. 
At smaller scales (<~ 7°), temperature variations at high 
latitudes (outside large molecular complexes) indicate a 
variation of the ISRF smaller than ~ 30%. 

It is also possible to build a map of large grain tem- 
perature with the FIRAS data by fitting each spectrum 
with a modified Planck curve. The signal to noise ratio 
of individual spectra is low and temperature uncertain- 
ties are typically two times higher than the DIRBE ones. 
However, differences of temperatures obtained from these 
two independent datasets are within the uncertainties. 
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Fig. 2. Temperature of the large grains (see Fig. 1) as a func- 
tion of the Galactic latitude. Pixels containing known com- 
plexes are identified using different symbols: Taurus, Auriga 
and Perseus (o), Orion (x), Large Magellanic Cloud (A), 
Ophiucus(*) 

The ISM is known to contain structures over a very 
wide range of scales with inhomogeneous physical con- 
ditions. Temperatures presented in Fig. 1 correspond to 
values averaged along the lines of sight and within the 
7° FIRAS beam. For pixels with a high signal to noise 
ratio (especially at |6| < 30°), the DIRBE data without 
any convolution allow the visualisation of spatial varia- 
tions of the temperatures at a resolution of 40' (Fig. 3). 
For the Taurus, Auriga and Perseus (TAP) region, the 
12 CO ( J = — 1) integrated emission map obtained by 
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Fig. 1. Temperature map of the large grain dust component computed from the DIRBE maps at 240 and 140 /jm convolved 
with the FIRAS PSF (in Galactic coordinates) . 



Ungerechts & Thaddeus (1987) at a resolution of 0.5° is 
compared in Figs. 4a-b to the temperature map at the 
DIRBE resolution (projected upon the 12 CO grid). We see 
that, within this field, there are significant temperature 
variations (Fig 4b) . The regions of cold temperatures spa- 
tially correlate with the 12 CO filaments. The highest tem- 
perature in Perseus (a — 3hA0ml0.6s, S — 31d58m7.9s) 
corresponds to the IC 348 HII region (also known as the 
Per OB2 cloud). 



4. Cold emission maps 



4-1. Maps of FIR excess with respect to the 60 /im emis- 
sion 

Studies based on IRAS images have revealed, for dense 
clouds, sharp increases of I„(100) without any counter- 
part of I„(60) (Laureijs et al. 1991 and Abergel et al. 
1994). It is illustrated for the TAP region Figs. 4c-d. Im- 
age of the component of the dust emission in excess at 
100 yum (Fig. 4e) is computed using J l/ (100) ea:C es S =/ ; y(100)- 
I„(60)xi?(100,60), where R(100,60) is the 7 I/ (100)/7 I/ (60) 
ratio assumed constant for the neighbouring emission 
("cirrus" component). The emission in excess has been 
found to be correlated with the emission of dense gas in 
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13 CO. In this section, we extend the IRAS excess com- 
ponent analysis to longer wavclcnghts using the DIRBE 
data. 

Let I D (\) (A = 100, 140, 240 fjaa) be the DIRBE emis- 
sion at 100, 140 and 240 /im respectively, and R(A,60) the 
I v {\) / I v {<oQi) flux ratio observed in "cirrus" clouds. The 
excess maps are computed at each wavelength according to 
the relationship: ir>(A) eX cess = IdW — R(\ 60) * 7/5(60). 

The 7d(A) vs Id (60) correlation diagrams (used to 
derive i?(A,60)) present a high dispersion due to real 
color variations across the sky (due to temperature ef- 
fects and/or abundance variations). However, the signal 
to noise ratio and the angular resolution of DIRBE data 



do not allow a proper determination of the spatial varia- 
tions of R(X, 60) over the whole sky. Therefore, we have 
decided to work with averaged values estimated by eye- 
fitting the Id (A) vs Id (60) correlation diagrams for all 
pixels at \b\ > 20° and |/3| > 20°, outside known molecu- 
lar complexes and the Magellanic Clouds. These values are 
equal to 4 ±0.7, 6±1 and 3.2 ±0.5 at 240, 140 and 100 ^m 
respectively. The error bars are due to the dispersion of 
the points. 

Using these three color ratios, we can also compute 
R(140,100) and R(240,100) (see Table 1). They are re- 
markably close to the slope of the Id (A) vs Id (100) cor- 
relation diagrams in Dwek et al. (1997). The discrepancy 
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Fig. 4. Taurus, Perseus and Auriga region: (a) CO J=0-1 integrated emission from Ungerechts & Thaddeus (1987) (b) large 
grain temperature map (the same as Fig. 3 projected on the CO grid) (c) 60 /im and (d) 100 /im emission. Cold component of 
the dust emission at 100 (e) and 240 /im (f) corresponds to bright features. They are due to a sharp increase of 7d(100) and 
Id (240) without any counterpart of Id (60). These features are very well correlated with the CO emission, confirming that, in 
this region, molecular emitting regions coincide with cold regions (Abergel et al. 1994) 
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observed for R(100,60) is due to the fact that we do not 
study the same part of the sky (|6| > 30° and \/3\ > 40° for 
Dwek et al.). With the same data points as ours and the 
same data preparation as Dwek, a value of 3.4 is obtained 
(Arendt, private communication), which is very close to 
our value of 3.2. Our color ratios differ significantly from 
those of Boulanger et al. (1996), especially R(240,100). We 
have checked that these differences cannot be accounted 
for by the fact that we do not study the same part of the 
sky. The color ratios depend on which map is used as in- 
dependent variable for the correlation (HI for Boulanger 
et al., 100 /im for Dwek et al. and 60 /im in the present 
study). 

We have used our estimates of R(X, 60) (given in Table 
1) to compute the excess component maps at 100, 140 
and 240 /im. These maps still contain a zodiacal residual 
emission (mainly coming from the 60 /im map) visible on 
large scale. It is difficult to produce all sky maps so we 
decide to concentrate our study on local features. Low 
frequency structures in excess maps are removed using a 
21°x21° median filter. The filtering also corrects for the 
large scale variations of the color ratios i?(A, 60) and for 
the isotropic contributions (as the CFIBR). 

4-2. Temperature of the excess component in the Taurus, 
Auriga and Perseus region 

The excess maps at 100 and 240 /im in the TAP region 
are presented Figs. 4 e-f. We have chosen 8 lines of sight 
to determine the temperature of the dust associated with 
the filaments in excess (Table 2, Ti). It is clear that these 
filaments correspond to cold regions with temperatures 
around 13 K. 

To test the effect of spatial variations from cloud 
to cloud in the color ratios R(A,60), we have derived 
the temperatures using excess maps computed with lo- 
cal color ratios derived in the vincinity of the TAP region: 
R(240,60)=8, R(140,60)=10 and R(100,60)=5. The new 
temperatures we obtain are very close to the ones deter- 
mined with the averaged color ratios (Table 2, Ti and T2). 
We see that variations in the color ratios do not affect the 
determination of the temperatures by more than 0.2 K. 

Finally, to check whether the temperature determina- 
tion may be biased by the 60 /jm removal, we compute the 
temperature using an "ON-OFF" method: a background 
emission is estimated for each pixel of the filaments at a 
distance of typically 3° (OFF position). This background 
is locally removed from the emission. Temperatures de- 
rived with this method are very close to the ones derived 
with the excess maps (Table 2). Therefore, we conclude 
that all bright filaments in the excess maps of the TAP 
region definitively correspond to cold filaments. We will 
see in Sect. 5 that we can generalise this result for the 
whole sky at \b\ > 10°. So, in the following, the excess 
maps are called "cold emission maps" . 



4.3. Correlation of the cold component with molecular gas 

The whole sky map of the cold component at 240 /im is 
presented in Fig. 5. We exclude in our analysis the regions 
at \b\ < 10° because of high spatial variations of R(X, 60). 
All bright features at |6| > 10° are associated with large 
molecular complexes: Polaris, Camelopardalis and Ursa 

Major (6 - 30°, I - 140°), Chamaeleon (6 16°, 

I ~ 300°), part of Orion (b 15°, I - 210°) and Taurus 

(b ~ —20°, I ~ 180°) as illustrated in more details in Figs. 
4e-f. The network of bright filaments at negative latitudes 
around b ~ —35° and I ~ 90° is associated with molecu- 
lar emission near the OB association I Lacerta. However, 
not all molecular complexes are visible on Fig. 5. For ex- 
ample, most of the p Ophiuchus complex is not detected, 
since it contains several B stars and a large number of 
embedded young stellar objects. We can also compare the 
240 /im cold map to the map of the residuals of the IR/HI 
correlation (Boulanger ct al. 1995). All structures in the 
cold maps appear as positive excess in the IR/HI residu- 
als, but the reverse is not true. Locally heated molecular 
clouds and ionised clouds are not present in Fig. 5 (as the 
one around the nearby Bl star Spica located at I = 316°, 
b = 51°, Zagury et al., 1998). 




5 10 15 20 25 30 
N HI * 10 20 (crrT 2 ) 



Fig. 6. Correlation diagram between the DIRBE 240 fim and 
the HI emissions convolved in the FIRAS beam. The data 
points for which the cold component is detectable at 240 um are 
surrounded. The continuous line, from Boulanger et al. (1996), 
represents the fit to data at N ff /<4.5 10 20 H cm -2 

Boulanger ct al. (1996) have correlated the FIR emis- 
sion from dust as measured by COBE at high \b\ with 
the 21-cm emission from gas as measured by the Lei- 
den/Dwingeloo survey of HI (Hartmann & Burton, 1997). 
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Fig. 5. Cold component of the dust emission at 240 /xm, 7d(240) coH = J D (240) - fl(240, 60) * J D (60), where i?(240, 60) is 
obtained by computing the slope of the Id(240) vs 7d(60) correlation diagrams for all pixels at |6| > 20° and \f3\ > 20°, outside 
large molecular complexes and Magellanic Clouds (see the text). We have removed the data at |&| < 10° because of high spatial 
variations of R/240. 601. 
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Table 1. FIR and submm color ratios for the whole sky at high Galactic latitude 



Publications 


R(100,60) 


R(140,60) 


R(240,60) 


R(140,100) 


R(240,100) 


This paper (1) 


3.2 ±0.5 


6± 1 


4 ±0.7 


1.9 ±0.6 


1.25 ±0.4 


Dwek et al. (1997) (2) 


6.25 






1.93 


1.28 


Boulanger et al. (1996) (3) 








2.36 


2.02 


(1) Ratios are estimated via the 60 /im emission (for \b 


> 20° and \(3\ > 20°) 



(2) Ratios are estimated via the 100 /im emission (for |6| > 30° and |/3| > 40°) 

(3) Ratios are estimated via the HI emission (for \b\ > 20°) 



Table 2. Temperatures of the bright filaments in the TAP region determined using (1) the excess maps derived with the 
averaged color ratios (Table 1 and Sect. 4.1), (2) the excess maps derived with local color ratios computed in the TAP region 
and (3) an "ON-OFF" method. 



RA (1950) 


DEC (1950) 


Ti (1) 


T 2 (2) 


Ta (3) 


04h03m33.64s 


40d24m36.8s 


14.6±0.3 


14.6±0.3 


15.3±0.6 


04 13 11.65 


28 43 12.6 


12.9±0.2 


13.0±0.2 


13.2±0.2 


04 16 15.12 


38 10 35.1 


14.0±0.2 


13.9±0.2 


14.3±0.2 


04 20 20.24 


37 24 39.3 


13.5±0.2 


13.3±0.3 


14.9±0.2 


04 23 41.15 


27 20 1.8 


13.9±0.3 


13.9±0.3 


14.9±0.3 


04 27 6.25 


24 40 18.0 


12.6±0.3 


12.5±0.3 


13.7±0.4 


04 27 18.82 


24 22 7.6 


12.5±0.2 


12.4±0.2 


13.3±0.3 


04 35 9.36 


26 14 4.6 


13.2±0.3 


13.0±0.2 


13.5±0.2 



At all wavelenghts (A > 100 /im), the correlation is tight 
for Nff/<5 10 20 H crrr 2 . For higher N ff/ , the data points 
depart from the low emission correlation. They interpret 
the change in the slope in the FIR-HI correlation as an 
increasing contribution of molecular gas for N HI larger 
than 5 10 20 H cm" 2 . The 240 /im-HI correlation at the 
FIRAS resolution is shown in Fig. 6 for |6| > 10°, using 
the same HI dataset. We have surrounded all data points 
for which the cold component is detectable (see Sect. 5.1) 
at 240 /im. Most of these " cold pixels" depart from the 
low correlation emission (continuous line of Fig. 6 from 
Boulanger et al., 1996) because of the increasing contri- 
bution of cold molecular regions. The other points which 
also depart from the low HI column density correlation 
correspond to emissions coming from (1) HII regions, (2) 
warm molecular clouds and (3) cold molecular clouds too 
small to be detected as such in the FIRAS beam. The min- 
imum HI column density for the "cold pixels" at 240 /xm is 
2.5 10 20 H cm~ 2 . We consider this value as a threshold for 
the detection of cold emission in DIRBE data convolved 
in the FIRAS beam. 



5. FIRAS spectra of the warm and cold compo- 
nents 

In this section, we combine the DIRBE and FIRAS data 
to extend the study of the Galactic emission (especially 
the separation between the warm and cold components) 



to longer wavelenghts. Thus, we work at the FIRAS reso- 
lution of 7°. 

5.1. Temperature and optical depth of the cold and warm 
components 

We convolve the cold DIRBE emission maps with the 
FIRAS PSF (I F (2A0) cold , I F (U0) cold and I F (100) coM ). 
These maps contain the emission of the cold dust com- 
ponent in the FIRAS beam. Even with the low FIRAS 
resolution, pixels for which the cold component is signif- 
icant can be identified. To select these pixels, we com- 
pute the root mean square (a\) for each cold map by fit- 
ting the brightness histograms by a gaussian curve. We 
find (7240=0.41 MJy/sr, eri 40 =0.54 MJy/sr and eri O o=0.22 
MJy/sr. The dominant sources of noise are spatial varia- 
tions in R(X, 60). We restrict our analysis to |6| > 10° and 
separate the data in two sets with and without significant 
cold emission. The pixels with significant cold emission 
are defined with: IfWcoI(1>3o-\ at 100, 140 and 240 /im. 
Mean cold emissions are equal to 1.0, 3.5 and 3.9 MJy/sr 
at 100, 140 and 240 /im respectively. 

The spectra of pixels with no significant cold emission 
(61% of the sky) are fitted by a single temperature Planck 
curve with a v 2 emissivity law. The mean temperature for 
this set of pixels is equal to 17.5 K with a dispersion of 2.5 
K. The high dispersion comes from low signal to noise ra- 
tio of individual FIRAS spectra at high Galactic latitude. 
Therefore, we confirm that the averaged temperature of 
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the interstellar dust in the atomic medium is rather uni- 
form and of about 17.5 K (Desert et al. 1990, Boulanger 
et al. 1996). 

For pixels with significant cold component (206 pixels 
which represent 3.4% of the sky), we measure the opti- 
cal depth and temperature of the warm and cold compo- 
nents. Let /f (A) be the intensity of the DIRBE emissions 
at the FIRAS resolution, IpWcirrus the intensity of the 
" cirrus" contribution (which represents the HI clouds and 
warm molecular envelops) and If{^)coI<i the intensity of 
the cold dust contribution. We compute If W cirrus using: 

iFWcirrus = -Tf(A) - IfWcoU 

The temperature and optical depth of the "cirrus" 
component (T c i rrus and T c i rrus ) are computed using 
If (240) c i rrw5 //F(140) C i rrus and assuming a v 2 emissivity 
law. Then, cold component spectra are obtained by remov- 
ing the "cirrus" modified Planck curves, defined by T C i rrus 
and T c i rrusi from the FIRAS spectra. Each of these cold 
component spectra is fitted by a modified Planck curve 
with a v 2 emissivity law (the C + line at 158 ^m is removed 
before doing the fits). We have checked that a free a does 
not significantly decrease the x 2 value of the cold com- 
ponent fits. As an illustration, fits of the cold component 
spectrum at I = -305.2° and b = -15.75° with (a) T=19 
K, t = 1.610~ 7 , a=lA and (b) T=15.25 K, r = 3.910~ 8 , 
a~2 give the same x 2 value of 1.26. We see that the poor 
signal to noise ratio does not allow to constrain a. More- 
over, the signal to noise ratio of individual cold spectra is 
too low for about 64% of them to even determine properly 
the temperature and optical depth assuming a=2. 
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Fig. 7. Histograms of the temperatures for the "cirrus" com- 
ponent (continuous line) and the cold component (dashed line), 
for pixels containing significant cold emission (see the text) 



The signal to noise ratio is larger with DIRBE than 
with FIRAS data. Thus, we use the DIRBE 100, 140 and 
240 /im cold maps to compute the temperature T co id and 
optical depth t co u of the cold component. We assume a 
v 2 emissivity law since the emissivity index (a) cannot be 
properly constrained. This v 2 emissivity law corresponds 
to standard interstellar dust grains (Draine & Lee, 1984). 
The cold component has optical depths (normalised to 1 
cm -1 ) distributed around 5.3 10~ 8 with a dispersion of 
4.1 10~ 8 and temperatures distributed around 15 K with 
a dispersion of 0.8 K (Fig. 7). The coldest temperatures 
(~13 K) are located in the Taurus region. Then, we derive 
spectra of the dust outside cold regions (" cirrus" compo- 
nent) by subtracting the cold modified Planck curve de- 
fined by T co id and t co u from the FIRAS spectra. Finally, 
each " cirrus" spectrum is individually fitted in order to de- 
termine T c irrus and Tcirrus- The "cirrus" component has 
optical depths distributed around 9.8 10 -8 with a disper- 
sion of 6.1 10 -8 and temperatures distributed around 17.8 
K with a dispersion of 1.2 K (Fig. 7). These temperatures 
are significantly different from the cold temperatures de- 
termined for the same set of pixels. An example of the 
decomposition in two components is shown in Fig. 8 for 
one pixel located in the Chamaeleon complex (I = 305.2° 
and 6= -15.75°). 
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Fig. 8. Example of the decomposition into the "cold" and 
"cirrus" components for one selected pixel (I = —305.2° and 
b = —15.75°) located in the Chamaeleon cloud. The pa- 
rameters of the modified Planck curves are: T co id~14:.2 K, 
t coH =6.9 10~ 7 (dotted line) and T cirru3 =19.2 K, r cirrus =5.5 
10~ s (dashed line). The solid line corresponds to the synthetic 
total emission. 
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5.2. Residual emission 

The presence of very cold dust (T=4-7 K), not detectable 
using the DIRBE data at 140 and 240 jUm, may be re- 
vealed on the residual FIRAS emission. We compute two 
residual spectra: one for regions with detected cold emis- 
sion, by removing the cirrus and cold contributions from 
the FIRAS spectra, and one for regions with no detected 
cold emission, by removing only the cirrus contribution 
from the FIRAS spectra. In Fig. 9, we present the mean 
emission spectra of the warm and cold components to- 
gether with the mean residual spectrum for regions with 
detected cold emission. The mean residue is negative in the 
range 280 — 390 /um (2% of the total mean emission) and 
positive above 400 /xm (~15% of the total mean emission 
at 600 /an). We have checked that this pattern persists if 
we fit the spectra with a fixed emissivity index different 
than 2. 




1000 



A (/xm) 

Fig. 9. Mean emission spectra of the warm (dashed line) and 
cold (dotted line) component for regions with detected cold 
emission. The continuous line shows the mean residual spec- 
trum (multiplied by 5), also for regions with detected cold emis- 
sion 

The residual to the one temperature fit of spectra with- 
out cold emission shows a similar pattern, comparable 
to the lcr positive submm excess of Dwek et al. (1997). 
This excess is typically 25 times smaller than the one de- 
tected in Reach et al. (1995). It could be due to very cold 
dust but we have checked that a temperature distribu- 
tion of the Galactic components qualitatively reproduces 
the spectral shape of the residues. We have seen (in Sect. 
5.1 and Fig. 7) that the temperature of the "cirrus" is 
~ 17.5 K while the temperature of the cold regions can be 



as low as 13 K. Obviously, inside the FIRAS pixels with 
detected cold emission, there is a range of temperature 
at least from T C i rrus to T co u- Moreover, the cold regions 
appear to coincide with molecular clouds which have a 
typical size smaller than the FIRAS and DIRBE beams. 
Therefore, pixels with no detected cold emission can also 
contain small size molecular clouds inducing a range of 
temperatures. Our residual emission does not contain any 
detectable very cold dust. 

6. Discussion 

6.1. Warm and cold components 

We have used the spectral information of DIRBE and FI- 
RAS to separate the different emissions of the ISM. The 
predominant emission of the whole sky in the FIR comes 
from "cirrus" clouds (as defined earlier: HI clouds and 
warm molecular envelopes). The large grain temperature 
in "cirrus" is relatively uniform and around 17.5 K (Sect. 
5.1). We do not see any systematic variation of the tem- 
peratures with the Galactic latitude which suggests a rel- 
atively constant ISRF on very large angular scales. For 
a v emissivity law, the dust temperature scales as the 
radiation field intensity to the power 1/6 and the spatial 
variations of temperature can be converted to variation in 
the radiation field intensity smaller than ~ 30%. The cold 
component of the dust emission with a temperature ^15 
K is spatially correlated with large molecular complexes 
with low star forming activity. The coldest temperature 
is <~ 13 K. Since these temperatures are averaged inside 
the FIRAS beam, the physical value of the dust temper- 
ature in molecular regions can obviously be lower than 
13 K. Such cold temperatures have also been found by 
the balloon-borne experiment PRONAOS in several dense 
cores in star forming regions with an angular resolution of 
2-3.5' (Ristorcelli et al., 1996, 1998, Serra et al., 1997). 

An important question is to know whether the drop in 
the temperature is only due to the attenuation of the radi- 
ation field in the molecular clouds or whether it also results 
from a change in the large grain properties. If the small 
particles disappear by sticking onto large grains within 
dense gas, it is natural to speculate that the emissivities 
of the large grains are different. The large grains could be- 
come fluffy and this will affect the long wavelength emis- 
sivity (e.g Bazell and Dwek, 1990). The IRAS emission at 
12 and 25 /im is observed to drop where the cold emission 
is present, and this has been interpreted as due to a low 
abundance of PAHs and VSGs (Bernard et al., 1992). We 
find here that low abundances are correlated with the drop 
in the temperatures of large grains. We derive an averaged 
value of R( 140,240) in cirrus clouds and in cold regions 
using the temperature histograms (Fig. 7). In Table 3 we 
compare these color ratios with those determined in the 
Bernard et al. (1992) model which computes with a radia- 
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tive transfer code the FIR emission of a non-homogeneous, 
spherically symmetric cloud with a standard dust compo- 
sition. The large error bars of our estimates of R(140,240), 
due to the dispersion of the temperature from clouds to 
clouds, do not allow to conclude if the drop in the tem- 
perature is only due to the attenuation of the radiation 
field or if it also results from a change in the large grains 
properties. It is necessary to study in detail individual 
molecular clouds which is beyond the scope of this paper 
which presents a statistical analysis on the whole sky. 



Table 3. Values of the R(140,240) ratio derived in cirrus clouds 
and in cold regions with DIRBE data and computed with the 
model developped by Bernard et al. (1992) for the case /3=1, 
central density n5j=10 4 cm' 3 and A„=l, 4, 6, 20. The value 
obtained by Desert et al. (1990) in cirrus clouds (A„=0) is also 
reported. 



,Md) 





R(140,240) 


Cirrus clouds 


1.3 ±0.4 


Cold regions 


0.9 ±0.4 


A„=0 


1.26 


A„=l 


1.1 


A„=4 


0.77 


A.^6 


0.64 


A„=20 


0.46 




Fig. 10. (a) Mean FIRAS spectrum in the region with 
30°<b<60° and 180°<1<270°. (b) Residuals after removing 
a single component modified black body with a=2. (c) Resid- 
uals after removing a single component modified black body 
with a=2 from the mean spectrum (a) corrected for the refer- 
ence zero emission spectrum (see Sect. 2 for more details). We 
see that, at high latitude, more than 80% of the submm excess 
(b) is due to the reference zero emission spectrum. 
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Fig. 11. (a) Synthetic spectrum (sum of 2 modified Planck 
curves with a=2, T coH =15 K, r coid =5.3 10~ 8 and T clrrus = 17.8 
K, r C irrus=9.8 10~ 8 ). (b) Residuals after removing a single 
component modified black body with a— 2 from the synthetic 
spectrum, (c) Residuals after removing a single component 
modified black body with a— 2 from the synthetic spectrum 
added of the reference zero emission spectrum. The residual 
emission contains the same submm excess as the one shown 
in the Fig. 3 of Reach et al. (1995) (d) Residuals after remov- 
ing two components modified black body with a=2 from the 
synthetic spectrum with the reference zero emission spectrum 
added (Ti=17.1 K, n=1.5 10~ 7 , T 2 =8.4 K, r 2 =8.15 10~ 8 ). 
The detection of the very cold component at 8.4 K is an arte- 
fact due to the combination of the cold one and the reference 
zero emission spectrum. The fluctuations in (b,c,d) are due to 
the standard deviations taken from the data (the same as in 
Reach et al.). 



6.2. Very cold dust in the Galaxy ? 

The problem of the presence of a very cold component 
in the Galaxy (T=4-7 K) is very important since a pos- 
itive detection could possibly imply the existence of an 
unknown population of Galactic grains. Moreover, a very 
cold dust component at high Galactic latitude would pro- 
duce a submm Galactic foreground extremely confusing 
for the future experiments dedicated to study the cosmo- 
logical backgrounds. In this section, we compare our re- 
sults with the previous studies of the submm emission of 
dust from balloon-borne experiments (Fischer et al., 1995 
and Masi et al., 1995) and FIRAS data (Wright et al., 
1991 and Reach et al., 1995). 

Fischer et al. and Masi et al. have detected a submm 
excess (with respect to a v 2 modified Planck curve fit) 
around the star ji Pegasi and in the Aries and Taurus 
region respectively. Their results are strongly constrained 
by the absolute values of the IRAS 100 /im emission. How- 
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ever, these values are uncertain, since the 100 /xm calibra- 
tion of IRAS is known to be different from that of COBE, 
and depend on the angular scale of sources. COBE/IRAS 
brightness ratio of 0.7 applies to very large scale emis- 
sion (Sodroski et al., 1994) but this factor is close to 1 
for angular scales smaller than 1°. Moreover, we detect a 
significant cold component with a temperature ~ 15.5 K 
in the FIRAS maps in both regions but no very cold dust 
emission. We can conclude that the regions mapped by 
those two balloon borne experiments arc not representa- 
tive of the diffuse interstellar medium and contain submm 
excess due to cold molecular clouds. 

Reach et al. (1995) have fitted FIRAS spectra, aver- 
aged in several bins over the whole sky and in the Galac- 
tic plane, with two components and found that the optical 
depths of the very cold component (4-7 K) and warm com- 
ponent (16-21 K) are correlated. This result strongly sup- 
ports the idea that the submm excess due to the very cold 
component has a Galactic origin. However, the Reach et 
al. analysis is based on total power FIRAS spectra which 
could include a non-Galactic emission component. Puget 
et al. (1996) have found, in the FIRAS data after the re- 
moval of the HI correlated emission, a positive residual 
emission. They argued that most of this emission comes 
from an isotropic component. At high latitudes, the pres- 
ence of such an isotropic component can explain the detec- 
tion of very cold emission with a two component fit (Fig. 
10, to be compared to Fig. 3b of Reach et al.). 

At intermediate latitudes (10° < \b\ < 30°), the 
isotropic emission becomes negligible and therefore cannot 
explain alone the very cold component found by Reach et 
al. However, we have seen that most of the cold emission 
(~15 K) is detected using DIRBE data in this part of the 
sky (Sect. 4 and Fig. 5). To check whether the submm 
excess of Reach et al. is due to this cold emission, we have 
built a synthetic spectrum (Fig. 11a) made of two modi- 
fied Planck curves with temperatures and optical depths 
equal to the averaged values found for pixels with detected 
cold emission (T co;<i =15 K, t coM =5.3 10 -8 , T cirT . us =17.8 
K and T C i rrus — 9.8 10~ 8 , see Sect. 5.1). The residual emis- 
sion after removing a single component modified black 
body shows a submm excess comparable to that found in 
Reach et al. (Fig. lib). The addition of the isotropic com- 
ponent to the synthetic spectrum increases the amplitude 
of this excess (Fig. 11 c). The fit of the synthetic spectrum 
with two modified Planck curves produces a component at 
8.4 K (Fig. lid), which is an artefact due to the combi- 
nation of the cold component at 15 K and the reference 
zero level emission spectrum. Therefore, we conclude that, 
out of the Galactic plane, the detection of very cold dust 
is due to the combination of the cold and isotropic com- 
ponents. The optical depth of the very cold component 
found by Reach is correlated to the optical depth of the 
warm component (Fig. 6 of Reach et al., 1995) since the 
cold molecular emission (^15 K) globally correlates with 
the warm emission (cirrus clouds). In the Galactic plane, 
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the detection of two correlated components is due to the 
high diversity of physical conditions (from HII regions to 
cold molecular clouds). 
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Fig. 12. Mean FIRAS spectrum in regions with significant 
cold emission (continuous line) is shown together with its de- 
composition into the warm (dashed line) , cold (dashed dotted 
line) and reference zero emission (dotted line) components. 
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Fig. 13. Mean FIRAS spectrum in "cirrus" regions is shown 
together with its decomposition into the warm (dashed line) 
and reference zero emission (dotted line) components. 
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7. Conclusion 

We have combined the DIRBE and FIRAS data to anal- 
yse the distribution of dust temperatures in the nearby 
interstellar medium seen at \b\ > 10°. The main results of 
this work follow: 

(1) The sky- averaged temperature of the dust in the whole 
Galaxy shows variations between lines of sight containing 
star forming molecular regions, clouds forming low mass 
stars, and cirrus clouds. At high latitude, the spatial vari- 
ations of the temperatures can be converted to variations 
in the radiation field intensity smaller than ~ 30%. 

(2) Maps of the excess of the FIR/60 /im correlation trace 
the cold component of the dust emission. This cold compo- 
nent appears to coincide with large molecular complexes 
with low star forming activity (such as Taurus, Polaris, 
Chamaeleon...). The association between the cold compo- 
nent and molecular clouds is further demonstrated by the 
fact that all sky pixels with significant cold emission have 
an excess IR emission with respect to the high latitude 
IR/HI correlation. 

(3) The combination of DIRBE maps of the cold emission 
with FIRAS spectra have pointed out the significant dif- 
ference in the dust temperature for the diffuse and dense 
parts of the interstellar medium (T distributed around 
17.5 K and 15 K respectively with a v 2 emissivity law in 
both cases inside a beam of 7°). 

(4) The FIRAS spectral residues of our analysis do not 
indicate any detectable very cold dust (4-7 K) present in 
the Galaxy at the level found by Reach et al. (1995). 

(5) We conclude that the "detection" of such a very cold 
dust in our Galaxy is due to the presence of an isotropic 
background (Puget et al., 1996), together with the cold 
molecular clouds identified from DIRBE data. FIRAS 
spectra are very well decomposed by one isotropic com- 
ponent and one or two Galactic components (Fig. 12 and 
13). 
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